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Abstract 

The theory of three dimensional turbulent wall jets 
with an infinite submergence is presented, followed by a 
brief discussion on how the theory can be adapted to the 
Situation of a finite tailwater. the tailiwater conditions 
were experimentally found to consist of three states. 
Firstly,shallow tailwater in which the theory is invalid. 
Secondly,moderate tailwater in which the theory is valid for 
a relatively short streamwise distance. Finally,deep 
tailwater state in which the theory remains valid. In the 
streamwise direction,the three dimensional turbulent wal] 
jet with a finite tailwater was found to consist of three 
regions,firstly the potential core region,secondly the 
radial type decay region and thirdly an approximately 
"uniform" flow region,with a transition zone between the 
second and third regions. 

An experimental investigation utilising a flow 
visualisation technique provided evidence that,except for 
very low tailwaters,the expansion of the jet was linear and 
dependent on the tailwater depth. The change of jet 
expansion with the change in tailwater depth was found to be 
described by a unique curve. From experimental measurements 
it was found that both the vertical and transverse velocity 
profiles were similar. The decay of maximum centreline 
velocity was shown to be less for shallower tailwaters than 
for deep tailwaters and using the concept of the virtual 


origin, the decay of maximum velocity has been presented in a 
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convenient form. The change-in the length of the potential 
core and the surface velocity were found to be dependent on 
the tailwater depth. The length scales for the velocity 
profiles were studied and curves for their prediction have 
been developed. A preliminary investigation of the bed shear 
stress has shown that the centreline shear stress decay is 
reduced as the tailwater becomes shallower and that the 
transverse shear stress profiles are similar. Finally the 
scour hole produced by the bluff wall jet eroding a sand bed 
was found to be similar and that the presence of a scour 


hole seriously affected the jet expansion. 
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I. INTRODUCTION 


A. Turbulent Jets 


When a jet of water enters a large body of water or a 
jet of air enters into a large expanse of air ,surfaces of 
discontinuity or separation appear .The separation surfaces 
are very unstable and result in the formation of eddies that 
move in a random fashion,both along and across the jet.This 
motion causes a transfer of momentum between adjoining 
havyerswane Fests sim aderiaite regi Gneundtaiids. 2 continuous 
velocity distribution and is Known as a shear layer.For 
large Reynolds numbers,as is the case in most practical 
Situations,this shear layer is turbulent.The turbulent 
mixing gradually extends into both the jet and the 
surrounding fluid causing on one hand,after a short 
distance,a fully turbulent jet and on the other, an 
expansion of the jet. 

If a jet diffuses into the ambient fluid unaffected by 
any solid boundary, the jet is Known as a free jet.Commonly 
investigated free jets are the plane jet and the 
axisymmetric jet. A plane jet is one issuing from a nozzle 
that is(theoretically) infinitely long,while an axisymmetric 
Jet SsOner1SssSuings trom decirculan MoZzz ie, 

The wall jet is defined as a jet of fluid issuing 
tangentially to a boundary surrounded by stationary (or 
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plane turbulent wall jet issuing into the same stationary 
fluid of semi-infinite extent on a smooth boundary and has 
been fairly extensively researched(Rajaratnam, 1976). 

The third type of turbulent jet,either free or wall, is 
the three dimensional jet.A jet issuing from a nozzle that 
is neither axisymmetric nor rectangular!with a large aspect 
ha tiOnwchat is, the mauio. of ethe width ormatnesslolsto the 
height)is Known as a three dimensional jet.Three dimensional 
jets can be free or wall jets but the only type to be 
considered here in detail is the wall fon 

Three dimensional! wall jets can be either slender or 
blutf wall gets .lfva get with a tairmly Wearqe aspect ratio 
issues tangentially onto a smooth plate,the jet is found to 
consist of three distinct regions. Near the nozzle there is 
a region called the potential core,in which the maximum 
velocity along the central plane,normal to the bed, remains 
constant and equal to the initial velocity.However as the 
shear layer on the upper side and the boundary layer on the 
bed grow ,the potential core is reduced until it disappears 
completely .From this point on downstream the maximum 
velocity decays. Initially the velocity decay is affected by 
the characteristics of the nozzle and was termed by Viets 
and Sforza(1966) the characteristic decay region.After a 
certain distance downstream of the nozzle the shear layers 
from the sides meet and the velocity decay becomes 
independent of the nozzle characteristics .ihis region is 


referred to as the radial-type decay because the velocity 
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decays approximately with 1/x,where x is the axial distance 
from: the out letuwhiicheis simi lar tosthesradialewall 
jet(Glauert,1956).A wall jet with these three pronounced 
regions is called a slender jet. 

ror aenozzie wiihkameaspect ratioewhich 1s close to 
unity ,it has been found(Sforza and Herbst,1970;Newman et 
al.,1972; Pani,1972)that the velocity decay passes from the 
potential core straight into the radial type decay 


region.This type of jet is classified as a bluff wall jet. 


B. Practical Considerations 


It is a common engineering problem to have a situation 
where a turbulent jet enters an ambient fluid of limited 
depth. Such situations occur when a pipe,culvert or control 
gate structure discharges into a shallow pond, lake or 
reservoir, especially in the littoral regions. The prairie 
lake is a notable case where the depth remains very shallow 
for most,if not all,of the lake.The Athabasca-Peace Delta 
region in northern Alberta is a prime example,where the 
depth is approximately two metres for the whole delta 
region,extending over an area of over 2,500km* 

(Kellerhals, 1971). 

In these situations the outlet will normally discharge 
on or very close to,and tangentially along,the bed of the 
lake. This problem can possibly be modelled as a turbulent 


wall jet with a finite tailwater,for the outlets described 
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above the aspect ratio will be small and therefore a bluff 
wall jet.In some cases ,such as a gorge,a river entering a 
resermvourtor lake may also bec lassif ted! as aya ludmr wall 

jetyailthough! in most cases “the aspect “ratio wall probably 


put a river into the slender jet range. 
C. Existing Work 


Although not many investigations of three dimensiona| 
wall jets have been undertaken, the studies that have been 
carried out are quite comprehensive. Viets and 
Sforza(1966)and Sferza and Herbst(1970)carried out 
experimental investigations of both slender and bluff wall 
jets .Pani(1972) carried out both a theoretical prediction 
and experimental work on bluff wall jets.Pani(1972) also did 
some work on circular wall jets with swirl.Newman et 
al.(1972) conducted a series of experiments with a circular 
nozzle resting on top of a smooth wall,while Chandrasekhara 
Swamy and Bandyopadhyay (1975) researched a circular air jet 
passing over a plate placed along the central axis but a 
short distance from the nozzle .Summaries of the state of 
knowledge on bluff wall jets are given by Rajaratnam( 1976) 
in his book"Turbulent Jets" and by Launder and Rodi(1981). 

Some work concerning erosion due to turbulent wall jets 
has been recently carried out.Rajaratnam and Berry( i977) 
studied the erosion due to a circular turbulent wall jet for 


the case when the tailwater was much deeper and the channel 
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much wider than the jet.Rajaratnam and Diebel(1981) extended 
this work to the case of very low tailwater and variation of 
channel width, while Rajaratnam( 1981) and Rajaratnam and 
Macdougal1(1981) have investigated erosion caused by a plane 
turbulent water wall jet with both deep tailwater and 
tailwater of the same depth as the jet thickness. 

The only published work on a turbulent jet with limited 
tailwater depth appears to be a study by Rajaratnam and 
Subramanya(1968)investigating hydraulic jumps below abrupt 
symmetrical expansions.In that paper the downstream channe| 
was of limited width,unlike the channel of the present 


study: 


D. Outline of Proposed Work 


Due to the lack of Knowledge on the study of three 
dimensional turbulent wall jets with finite tailwater 
depth, the experiments carried out were of a very exploratory 
nature.In order to Keep the investigation simple and prevent 
confusion in the data due to different jet conditions, the 
experiments were limited to a study of bluff nozzles.Also, 
as will be discussed later, the downstream distances 
involved in the work are,for jet flows,relatively short,only 
to a distance of approximately thirty times the nozzle 
height. 

The experimental work was divided into three sections. 


Firstly a flow visualisation study of the change on the 
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expansion of the jet due to the change of tailwater depth 
was investigated, with tailwater to nozzle thickness ratios 
between zero and approximately fourteen(approaching infinite 
submergence condition), using five different aspect ratios. 

Secondly a series of velocity measurements, both 
vertical and transverse were taken,along with bed shear 
stress measurements .The results were analysed with respect 
to similarity of profiles, velocity and shear stress decay 
and surface velocity,along with any related length 
scales.Tailwater to nozzle thickness ratios varied from 
approximately one to about five.Two aspect ratios were used. 

Thirdly a smaller piece of work investigated the 
erosion of a fine sand bed due to slow moving jets.The 
centreline scour hole profile was measured and analysed for 
Similarity. The tailwater to outlet depth ratio was 
approximately one and only one aspect ratio was used. 

A short discussion on the theoretical analysis of a 
three dimensional turbulent wall jet and the difficulties in 
adapting it to the case of a finite tailwater is also 
presented 

The experimental investigations were carried out at the 
T. Blench Hydraulics Laboratory of the University of 
Alberta, Edmonton,Canada during the spring,summer and autumn 
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E. Concept of the “Uniform Flow Region” 


As discussed earlier a bluff wall jet issuing into a 
semi-infinite stationary body of the same fluid has two 
distinct regions,the potential core region and the 
radial-type decay region.However it was found,during the 
course of the investigation ,that this was not the case for 
the situation of a finite tailwater and in order to promote 
a better understanding of the problem, a qualitative 
description of the jet will now be presented. 

Once the jet has entered the shallow body of fluid it 
is obvious that the jet will expand as discussed earlier.It 
is also easy to envision that due to the limited tailwater 
the jet will quickly fill the whole depth of taiiwater.Once 
the potential core region has disappeared ,the maximum 
velocity will begin to decay in the normal radial decay 
Fashion. However, as soon as the jet fills the complete 
depth the only stress decelerating the upper surface is the 
shear stress due to thelin this case)air-water 
interface,which is negligible compared to the flow 
deceleration and that due to the bed (or wall) shear stress. 
Therefore, the maximum velocity is decaying faster than the 
surface velocity and eventually a region of (approximately) 
uniform vertical flow is reached. 

In@orderstombass from therverticalmvelooity pronmite 
found in the radial type flow region to an(approximately) 


uniform vertical velocity profile a transition zone is 
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formed. 

As the tailwater depth decreases, then the jet expands 
more rapidly into the whole depth and therefore the 
"uniform" flow transition zone encroaches upon the radial 
type decay region , thereby shortening the potential core 
region.When the tai lwater depth to nozzle height ratio 
approaches unity the jet seems to expand abruptly into the 
whole depth and the radial type decay region no longer 
exists.The jet appears to pass from a shortened potential 
core into the uniform flow transition region. 

It is obvious that once in the the uniform flow region, 
the normal methods of analysing the similarity of velocity 
profiles will no longer be correct because the maximum 
velocity will occur at an infinite number of locations,even 
towards the end of the transition zone the choosing of the 
correct position of the maximum velociy would probably 
become difficult and this is one of the reasons why 
measurements were restricted to a fairly short downstream 
distance. 

No attempt was made to measure velocity profiles when 
the tailwater to nozzle ratio was less than unity and no 
attempt will be made to predict the flow patterns, these are 
areas which require further detailed investigation. 

Figures 1.1 and 1.2 give a schematic representation of 
the ideas proposed.Reference should also be made to the 
veriical velocity profiles obtained from Run1 and 
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Presented within this chapter is an outline of the 
theory of three dimensional wall jets issuing into an 
infinite half space of the same fluid. This is followed by a 
discussion on how the theory may be adapted for the 
condition of a finite tailwater depth and how under certain 
circumstances the assumptions made in the theory are 
invalid. Also presented is a short discussion on the use of 
dimensional analysis in this study. In this chapter, and the 
study as a whole, the cartesian system of co-ordinates is 
used. The x axis is along the bed and central plane of 
symmetry of the jet, the y direction is normal to the bed 
ano tne 2 axis 1S im the third or iunogenalvairectionm, "as 


Shown an tigure 2.115 


B. Equations of Motion 


The Reynolds equations for steady flow in the cartesian 


co-ordinates can be written as(Schlichting, 1968) ; 
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where x,y and z axis are defined in section A 


u,V,w=the turbulent mean velocity components 


x,y and z directions respectively 
ur wi?! 
in x,y, and z directions respectively 
p=mean piezometric pressure at the point of 
consideration 


o =mass density of the fluid and 


v =Kinematic viscosity of the fluid 


he 


»w =the fluctuating component of the velocities 
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C. Deep submergence Theory 


Presented in this section is a brief review of the 
theory of three dimensional! turbulent wall jets. The theory 
is not dealt with in great detail and only the salient 
points are noted: for a detailed discussion the reader 


should refer to Pani(1972) or Rajaratnam( 1976). 
Simplified Equations of Motion 


The following assumptions are made; 

(fi Beyond, an initial (reaches from: the nozzle y<<x and ez<<x 
aneiivees major portion ot thes tilow, Uo Vv andea ow: 

(7ii)The fluctuating velocity components are small compared 
to the mean values. 

(iv)Gradients in the y and z directions are much larger than 
in the corresponding gradients in the x direction. 

(v)The piezometric pressure gradients in all three 
directions are small. 

In estimating the order of magnitude of each term, the 
order of magnitudes for the dimensional forms of equations 
2.1 to 2.4 were assigned as 

U~1l ,x~1 ,y~a and Z~aA 
where Asis auquantTlyeofmorder much lessy Thansuniuhy ite is 
possible to icapny, outm@anyorder of magnitude analysis on the 
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the equations of motion are reduced to 
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Similarity Analysis 


If, when the ratio of the local velocity to the local 
maximum velocity is plotted against the ratio of the 
distance from the jet axis to some characteristic length, 
the velocity profiles collapse onto a single curve, the 
velocity profiles are termed "similar". In a wall jet the 
distance in the y direction is measured from the wall (or 
boundary}. The characteristic length has usually been eae 
as the distance to where the velocity becomes half that of 
the local maximum and a and 91s) Known as thes shalt 
length”. 

Glauert(1956) has shown theoretically that complete 
similarity of velocity profiles is not possible due to a 
thin layer existing near the wall where viscosity is 
important.However the experience of many investigators has 
shown that a single velocity scale and a single length scale 
correlates the data well. Therefore while it is not 


theoretically feasible, it is possible to have experimental 


at? fo. piten Sar sacte 41 
1 tM Tehtebé Bartels of er lespiev etal 
2 wt BUA fe ed many eonsd 


Girt eunklico asl tes ‘ginel ay edd 
\ yragon! Pet av 


1? brattiegant 26 Aol*oes'D y act? oF leone tate 

di je! -ei panes. ant | wneud “a 

a4 via ©! ss riedete eh! 2s. 7 
an stl) ef ncn. a! tre . D> . brey ion gem), Peered «it 

; a dtpnet - 


i ?aQr 3v i607 


e 


5 
isi int= beams! ove ca! 
iT? 
fiuew abet’ Wo oie 


ssmatad: VrisétivesAd a 


lapiignesds mathe ae (32er | Puehe 


eisai ?arnit yi 
® of su elatesog 104 2? 22! f honGivii(galsy Fe tant inti 


ai wi Raptety eee |) ew off seen ontietee vee nit 


D Sam arataptrseuni \iben Ng pSoriet rae aid ee 
| ipa Sael ate) 1 on ahead Signie we 


4 


Simplari tye ltssnout.ca bemnotedathat simian pron les are 
independent of the x co-ordinate. 


ASSUMING surat SROMm Samia ty 


Uu 
mo y 


where Um= turbulent mean velocity in the x direction at a 
distance y from the bed in the z=0 plane 

Unon ene local maximum velocity 

by =. New meal i eleng tne imme tnen yad weer ar 

f =a Tunctional relation 


amMdethat 1m the x=zeplane 


Ube (ees (2.8) 
Tene hy Pade 


where u=turbulent mean velocity in the x direction at a 
distance z from the axis 

bi=thesna its lengthen sthe Zz direc ion 

g, =a functional relation 
From experience with other jet flows it is possible to make 


the assumptions (Rajaratnam, 1976) 
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where g denotes a functional relationship. 


Ly it ees) hUmthereasstimed that 
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where p, q,and q, are unknown exponents. Then it can be 
shown that by manipulation of the above set of assumptions 


Eq 2.5 becomes 


Bei - ee Sioa) Bay a 
“Gh = g ie ; 1 95 
ele ee Ui Py 1 Um Py DS 19) bD 
b ou 1 
yom De eee : Ve 
Se i ak oR a SAL ea 
Mg yi 
Boe b ' 
Vette zy 
a tne ae pus eo oe nee a 8 


where, for convenience, (y/by) and (z/bz) have been dropped 
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also the prime on the f function denotes differentiation 
Witterespect to 1y/by) and on the g functions denotes 
Gifferentiation with respect to (z/b,). 

In turbulent jets teaby is generally small and so the 
iMiestoe term on the mightahand side Of “EG =Zaimecouldebe 
neglected, this “assumption will onlyerun anto dithiculties 
near the wall. Also for wall jets by/b, is less than unity, 
so that the second term on the right could also be dropped. 
For similarity, since the left side is independent of x, the 


right must also be independent of x. Therefore 
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Therefore 2 Bh. aaa DUttne value’ Of exponemiin pals sii || 


unknown. 


Momentum Integral Equation 


To evaluate the exponent p the momentum integral 


eGudtion can be utilised Fromocas 2 jo and Zao 


Z 2 
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The boundary conditions for deep submergence are, 


at y=0 u=v=w = 0 
u'= v'= w'= 0 
=) ee} = oe 
y U 0 yy 0 
ut=" 0 
and at 7a I) anise uo => 0 —_ 0 
7 4 
u'= 0 


Eq 2.16 is integrated over the infinite half-space bounded 


by y=0 plane and reduces to 
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is tne shear stress on the boundary 


whereé Tj 1S 
Rewriting the right hand side of Eq 2.17 as 


1 T 1 0 
- — dz = -—T x (——) dé (2218) 
0 | 0 0 Om es 
where x=distance measured from the nozzle in the x 
direction 


Tom =boundar y shear stress on the centreline and is 
in the x direction 
B= 7 x 


and assuming that 
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where, Chere alate! (GL empirical constants found from 
experimental data. Ug is the efflux velocity and A is the 
GrOss sect ionda |¥areal ome thennozzile. Eq 2e2siseonly a 
simplified version that does not fit the experimental data 
very well. 

Then it is possible that by substitution and 


integration to obtain the equation 
dz = [ mew ts a ] (2222) 
Ab ou-dydz Wg |Lt > eos LL. z, 


where My = eAUS =eF Flux momentum of the nozzle and Xg=the 
VaIUevOrex ate the Virtual sori. 

Based on his experimental results Pani(i972) was able to 
evalUateva je, and. Casand Simp iimy Faq 2. 22ePand i 1o7 2) then 
demonstrated that the momentum loss in a distance of one 
hundred times the nozzle thickness was around 11%. Therefore 
for a rough approximation the three dimensional wall jet 


momentum is conserved and Eq 2.22 reduces to 


2 
i fo dydz = Mp (2723) 
Exco 0 
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ed d (Mg) 
dx pu Cdydz Stee = 0 (224) 
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Therefore for a three dimensional wall jet 


and boa x (2.28) 
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Dimensional Considerations 


The use of dimensional analysis always involves 
engineering judgement and experience in order to choose the 
important variables.For a three dimensional wall jet with 
deep submergence the following results have been 
successfully vindicated by experimental observations 
Pamir on2.)4 

Generally, Reynolds numbers in jets are over a few 
Pecan and therefore it is reasonable to neglect the 
effect of viscosity. Therefore the maximum velocity along 
the centreline is a function of the initial momentum 
- flux(because of conservation of momentum), the fluid density 
and the distance X measured from a virtual origin(refer to 
chapter 4 for’ detailed discussion). 


= ( X 
en, O(M,,0,x) (2.29) 


Using Buckinghams Pi theorem it is possible to obtain that 


Uno a6 (2.30) 
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where C is an empirical constant 
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The centreline shear stress can be expressed as 


Bee = CHL pesky 2D) (2.32) 
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FOnue Given Reynolds number , tomo x2 
Because the snear stress is a function of the jet 
velocity, jt is reasonable to expect- that if the transverse 
velocity profiles are similar then the transverse bed shear 


stress, profides are also similar. 


no} ¥ '(2e) (2.34) 


and also that Boe the shear stress half length, x. 
D. Adaptation to Finite Tailwater 


The following paragraphs attempt to show how the theory 


of three dimensional wall jets in some cases remains to a 
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reasonable approximation valid, while under other conditions 
some of the assumptions made in the theory are no longer 


viable and the theory breaks down. 


Shallow Tailwater Condition 


From observations made of the flow, it became obvious 
that for tailwaters approaching the same depth as the nozzle 
height (and tailwaters less than the nozzle height) the 
surface of the water no longer remained flat. The waves and 
ripples formed on the surface were unsteady , sharp-crested 
and localised.Therefore, for a shallow tailwater it can no 
longer be assumed that the pressure gradients are 
negligible, nor. that any gradients in the x direction are 
smaller than those in the y and z planes. It is also 
possible that fluctuating velocity components may no longer 
be small compared to the mean velocity. These factors lead 
to the conclusion that the equations of motion can no longer 
be simplified as discussed above. Because the flow pattern 
has become so complex no attempt has been made to simplify 
the equations of motion and the theory for infinite 
submergence can no longer be accepted. 

The actual definition of a shallow tailwater is unclear 
and ,from observations ,is dependent on the Froude number fF, 


defined as 


mA . UG (2535) 
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where Vg Rilevere © height 

However, for the range of Froude numbers considered in 
this present study it appears that the surface disturbances 
become small for a tailwater one third greater than the 
nozzle height. It should also be considered that for a 
shallow tailwater the jet expansion in the transverse 
direction becomes very small(see chapter 4) and therefore 
the ratio of by/b, can no longer be assumed to be much less 


than unity, further invalidating the theory. 


“Uniform Flow Region" 


In chapter 1 the concept of an (almost) uniform flow 
region was discussed and it was stated that within this 
Pegion the similarity ot aventicaieve locityyprog? testwou la 
become undefined. It was also discussed that this region 
occurs in both shallow and moderate tailwaters. Therefore, 
at a certain distance downstream of the nozzle outlet, the 
similarity analysis used earlier can no longer be true and 
again the’ theory Tails. Again the definition ef a moderate 
tauiWalter 1SmUnG ledimancece deDenCel UlC nmr um OU lating ine 
present study a tailwater to nozzle height ratio of 
approximately four appears to be a reasonable upper limit. 
The distance downstream to where the (approximately) uniform 


flow region starts is undefined and 1s dependent on the 
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tailwater depth and probably the Froude number. It should 
also be noted that even for quite deeply submerged flow 
conditions, the jet expansion still reaches the the water 
surface. When the jet expansion reaches the surface some 
surraces disturbances dovroccur but sare notithought  toabe 


Sigmaificant. 


Moderate and Deep Tailwater Conditions 


Within thewremaininger lowes tualions#y imal ise) ror ea 
moderately deep tailwater, fairly close to the nozzle outlet 
and for the deep submergence condition it is possible to 
investigate the theory of three dimensional turbulent wall 
jets to discover whether it can still be used. 

The assumptions made in simplifying the equations of 
motion are still valid, as is the similarity analysis 
DEO Ginga hater Gm2 30.0) Ss ModiitedtsinGiul VeameOoria tant te 


tailwater it is reasonable to assume that 
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-but that the coeffecient of proportionality is different 
for each tailwater. That is, Eq 2.10 is true for any given 
tailwater depth. This, therefore, leads to the same 
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variation of the length scales are linear with the distance 
downstream from the nozzle but the equation of the line 
relating the two variables changes with the variation of the 
tailwater depth. 

Considering the momentum integral equation, this 
equation remains valid but the boundary condition for y= is 
no longer correct. At the water surface there is an 
air-water interface, on this interface it is reasonable to 
assume “no-slip" condition and within the bounds of accuracy 
used in the above work the shear stress on the interface can 
be neglected. Therefore the boundary condition for Eq 2.16 


becomes 


and the approximation of conservation of momentum is still 
correct. So again U,6%1/x tor amy given tailwater but, as 
before , the equation of the line is a function of the 


tailwater depth. 


Dimensional Considerations 


As stated in section C the use of dimensional analysis 
involves the use of considerable engineering judgement .One 
obvious approach is to use the Knowledge of three 


dimensional turbulent wall jets with infinite submergence as 
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a starting position and to make adjustments to account for a 
finite tailwater. 


For example the centreline velocity decay can be 
considered as 


(2.36) 
where y% =tailwater depth 
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theorem gives the dimensionless groups 


F M ; y 
pe }— ae (2.37) 
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Similarly the centreline shear stress could be expressed as 


0? Os 


Xo Vs Vas yee) (2.39) 
and, as shown earlier, 


for a given Reynolds number 


Tom! hi xX, 
but now the proportionality is a function of the tailwater. 
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The transverse shear stress profile can be expected to be 
similar but again the half width, b* , Vii I SCwnaaruhnculon OF 
iNew ls tance xsandethe Taniwater depth: ym 

It must also be recognised that the correctness of 
dimensional analysis is very sensitive to the variables 
Gn Ging! tyathoughtii to beminvelVved an themrUncrtOnm site, Tor 
instance, the wrong length scale is used then the 
dimensionless groups will be formed around the incorrect 


assumption. 
E. General Discussion 


As has been shown, the theory of three dimensiona | 
turbulent jets with a finite tailwater depth is not complete 
and in some areas non-existant. It is therefore necessary to 
complement theory by empirical relationships. This 
development of empirical eer Meets is attempted in the 
following chapters, using the results obtained from 
experiments performed on bluff jets with a finite tailwater 
depth. 

It must be recognised that experiments can only cover a 
limited range and this is true for the studies presented 
here. Also the development of empirical relationships 
involves a great amount of judgement, experience, trial and 
error and dimensional consideration. Therefore, the 
relationships and correlations presented in this study must 


be treated with caution and may need to be adjusted as more 
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data become available. Nor should they be applied to 
Situations which have different boundary conditions, for 
example when the jet is slender or when the jet discharges 
into a channel of restricted width. However the empirical 
results contained in this study do provide a far better 
estimate than trying to- apply infinite depth theory to 
finite tailwater cases and a better estimate than anything 


that, to the author’s Knowledge, is currently available. 


F. Summary 


An outline of the theory for infinite tailwater has 
been presented. The equations of motion were simplified 
using an order of magnitude analysis. Then using a 
similarity analysis and the momentum integral equation it 
was shown that the jet grows linearly and the centreline 
velocity decay is proportional to the inverse of the 
distance from a virtual origin. Using dimensional 
considerations it was shown that the centreline shear stress 
decay is proportional to the inverse of the square of the 
distance from a virtual origin for any given Reynolds 
number . 

It was then discussed how the theory became invalid for 
tailwater to nozzle height ratios of less than approximately 
13 and how the theory became invalid when the jet entered 
the (approximately) uniform flow region. Finally the theory 
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dependence on the tailwater, for flow situations in which 
the tailwater was deep, i.e. a tailwater to nozzle height 
ratio greater than approximately four or in moderate depth 


tailwater when restricted to relatively close to the nozzle. 


III. EXPERIMENTAL APPARATUS AND MEASUREMENT TECHNIQUES 


Atel mM troduction 


This chapter describes the apparatus used in the study 
and discusses the various measurement methods used. The 
limitations due to the apparatus and experimental procedures 
are also dealt with. The experiments were carried out in two 
flumes. The first flume, called the jet tank, was the 
smaller of the two but more versatile. The second and much 
larger flume, referred to (for the want of a better 
description) as the lake model, had a fixed tailwater 
condition and a free surface inlet (nozzle). The 
photographic study was carried out in both flumes, while 
velocity measurements were restricted to the jet tank and 


the scour investigation to the lake model. 


B. Jet tank 


The jet tank was of a rectangular shape, 3.55m long, 
1.11m wide and 1.22m deep, with plexiglass sides. A false 
bed was provided approximately at half depth. This bed was 
provided by a 5mm thick aluminum plate, 1.10m wide and 2.44m 
long, stiffened by two channel sections, running lengthwise, 
26cm from each edge. The plate was suspended at the upstream 
end by two, threaded, 5mm diameter steel rods, one in each 


corner. The rods passed through flanges attached to the top 
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of the tank ,thus enabling the upstream end of the plate to 
be leveled for each nozzle used. The downstream end of the 
plate was supported by hollow concrete construction blocks 
and a system of wedges, which allowed the plate to be set at 
the desired level. The test section of the bed was sprayed 
with white paint and a itcm grid was etched into the paint. 
running parallel and normal to the nozzle central axis. 

The inflow was fed, primarily from the laboratory sump 
but was supplemented by the city water supply when 
necessary, into a 150mm diameter standpipe. The standpipe 
could be varied in height to provide different constant heaa 
conditions. An overflow box collected the waste water and 
returned it to the sump. The standpipe turned through 90 
degrees into the horizontal plane and after a length of 76cm 
entered the nozzle. In the horizontal section of the pipe 
"hogs hair" was positioned to act as a baffle and reduce any 
fluctuations or secondary flow. The nozzles were well formed 
and the contraction took place over a 30.5cm length. A 
circular nozzle of 2.54cm diameter and a smail rectangular 
nozzle 0.64cm high by 1.58cm wide were turned from solid 
brass cylinders .A larger nozzle 2.5cm high and 4cm wide was 
shaped out of wood and sealed against water. 

The tailwater conditions were controlled by an 
adjustable standpipe 13.7cm in diameter. The outflow was 
collected and either returned to the sump or, if dye was 


being used, into the city sewer system. 
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Figure 3.1 gives a schematic diagram of the jet tank 


arrangement iandi plates Sad, 322, ses—ese4—snowsdetaids of 


the actual experimental setup. 


C. Lake Model 

The second flume was constructed from hol low 
construction blocks and was 15.25m long, 4.57m wide and 
0.61m deep. The inflow was provided from the laboratory sump 
and entered at the bottom of a head reservoir which was 
0.6im long, 0.405m wide and 0.6im deep. From the reservoir 
the inflow passed over a V notch weir and into a head canal 
0.9im long, 0.405m wide and 0.405m deep. The head in the 
reservoir was measured by a staff point gauge with an 
accuracy of 0.025mm and the discharge could be accurately 
calculated from the calibration curve of the weir 
(Hydraulics laboratory records). The inflow entered the lake 
model through a well formed constriction. The nozzle was 
50mm wide and the upper surface was unrestricted. The water 
level in the model was controlled by sharp-edged tailgates 
that were sealed at a height of 0.406m.The water returned to 
the laboratory sump, and a plugged 15cm diameter pipe passed 
underneath the tailgate, along the bed of the model, 
allowing the model to be drained if necessary. 

A 2.4m long by 2.4m wide sandbox of 15mm plywood was 
placed around the entrance to the model, the top of the 


sandbox was level with the bottom of the nozzle. The sandbox 
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Plate 3.1 - General view of jet tank arrangement 
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contained 0.1imm-mean diameter fine sand. For the 
experiments not involving scour investigations, a 3.66m long 
by 1.22m wide aluminum plate was placed over the sandbox and 
the entrance constriction was reshaped with plasticine. The 
downstream end of the plate was supported using hollow 
concrete construction blocks and wedges to ensure that the 
plate was level. 

Figure 3.2 shows a schematic representation of the lake 
model and Plates 3.5 and 3.6 show details of the actual 


arrangement. 


D. Photographic Procedures 


The photographic work was carried out in both flumes, 
and a Pentax MX 35mm camera was used througnout. Kodacolor 
Tl C135 print 11m was Used with an AS49100) settingsana a 
shutter speed of 1/30 second. 

For the jet tank series of experiments, the camera was 
supported by two horizontal angle bars running across the 
width of the tank. This enabled the camera to be set at a 
constant height of 1.5/7m above the bed for all the 
experimental runs. A food dye injection technique into the 
head pipe coloured the jet, allowing a photographic study of 
the jet expansion in) themitransVvemserdinect ion to besmade. 
Intially the dye was fed, under gravity, along 5mm diameter 
Dlastic UUbING tt mem SC alo OWe VehLin | Se OnOved 
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Plate 3.5 - Gemeral view of lake mode! 
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amount of dye to be injected under pressure. Intially it was 
hoped to take time averaged photographs with a slow shutter 
speed but due to the rapid discolouration of the tailwater 
by the food dye it was found necessary to increase the 
shutter speed. In order to obtain a true representation of 
the jet, at least three photographs were taken for each run 
and the values measured were averaged. The icm grid etched 
onto the bed allowed an accuracy of +2mm but there was 
considerable judgement needed in estimating the position of 
the nominal jet boundary. The jet exit velocity was measured 
wusitngea pitot-static ttubes which wilt be discussedmiater in 
detail and the tailwater was measured using a point gauge 
mounted to a traverse that was capable of moving in both 
horizontal directions. The point gauge could be moved 
vertically by a mechanical cog system, the gauge could be 
read to an accuracy of 0.025mm but play in the cogs reduced 
the accuracy of the reading to approximately 0.1mm. The lack 
of camera lense angle and the need for a clear grid, meant 
that the section of the tank photographed was quite short. 
For the lake model the camera was mounted on an angle 
beam, that in turn was bolted to a 3.5m high step ladder. 
The ladder was positioned over the head canal so that the 
camera was positioned approximately 3m above the sandbox 
test section. A scale was obtained by constructing a 10cm 
grid, of elasticated string with an aluminum frame, as shown 
inePhate S569 Thetorad waspphotographed on top of the 


sandbox and the metal plate before being removed for the 
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experimental runs. Because the camera was kept in the same 
posh onset ors alll Gehe mumset ther grad! photographs weouldebe 
used as a scale for all the photographs obtained from the 
study. Dye colouration was achieved by adding the food dye 
into the head canal. Again at least three photographs for 
each run were taken and the average value calculated. The 
mean outlet velocity was found from the discharge obtained 
from the V notch weir and the area of the nozzle. The depth 
of tailwater and the nozzle depth were measured using a 
point gauge, similar in construction to that described 
earlier for the jet tank and was mounted on a traverse 
capable of movement in both horizontal directions. The 
traverse could be positioned to +1.59mm in both directions. 
Plate 3.5 shows the camera mounting apparatus, while Plate 


3.6 illustrates the point gauge. 


E. Velocity Measurement 


Velocity readings from the jet tank were obtained using 
a comnercialby availablesPrandt|l=type "prtot static tube and 
commercially available pressure transducers. The pitot 
static tube had an external diameter of 1.59mm (1/16 in) and 
a hemispherical nose. The tube was clamped to the staff 
gauge on the traverse described above and is shown in Plate 
3.2. The tube was positioned manually and the etched grid 


was used as a position reference. 
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The pitot-static tube was connected to the pressure 
transoucer Using» 5mm diameter (plastic tTubingm Miegomessure 
transducer used for lower heads was a Celesco type P3900, 
whien had a full seale der lection of "#2. S4cm waters al though 
when tested it was found that the readings were stil] 
accurate to +8.0cm. For larger heads (greater than 2.54cm 
water) a Validyne model KP15, with a diaphragm full scale 
deflection of +0.5psi was used. This transducer had a full 
scale deflection of +35.0cm of water and, according to 
ManUtacturers specification, an accuracy of 075% anduan 
ability to withstand an overload of 200%. The system also 
used either a Pace CD25 transducer indicator or a Validyne 
CD15 carrier demodulator, both linked to a digital indicator 
that had two decimal places and was therefore capable (okg 
reading to 207@imV. Plate 3.3 ‘shows the pitot-static Lube, 
the Celesco transducer and the Pace CD25 transducer 
indicator and digital indicator. The apparatus was 
calibrated using calibration reservoirs, also shown in Plate 
3.3, the movable reservoir could be positioned to within 
+0,.025mm. The whole network system was filled with de-aired 
water that was coloured to facilitate the easy detection of 
air bubbles, which would invalidate any observations. A 
diagram of the system is shown in Figure 3.3 

It was usually necessary to take an average value for 
the velocity head due to the fluctuations in the flow. For 
large velocity heads, as was usually the case, this woula 


lead to negligible error but in areas of very slow flow an 
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error of 0.01mV could be quite significant(greater than 504) 
and in the “uniform flow region" of the jet, with small 
differences in velocity in the vertical plane, the 
Fluctuations of the readings could lead to significant 
error. This was another reason why. the investigation only 
studied the jet relatively close to the nozzle. The low heac 
transducer was calibrated to give a reading of 0.1mm of 
water for every 0.01imV, while the high head transducer gave 
a reading of 0.508mm of water for every 0.0imV shown on the 
digital indicator. It should also be noted that the 
pitot-static tube gives errors of less than 1% in velocity 
for misalignments to the direction of flow up to about 15 
degrees (Massey, 1975). 

As shown earlier, the velocity in the lake model was 
calculated by an indirect method. The discharge was found 
from the V notch weir and the nominal cross-sectional area 
of the nozzle was computed. Using the one dimensional! 
continuity equation the mean intial velocity of the jet 


could be calculated. 


F. Bed Shear Stress Measurement 


A very simple way to measure the boundary shear stress 
is to use a Preston tube(Preston,1954). A Preston tube is a 
total head tube resting on the boundary, for which the 
difference between the dynamic pressure on the tube and the 


static pressure on the boundary is correlated with the shear 
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stress. It has been shown (Rajaratnam, 1965; Rajaratnam and 
Muralidhar,1968) that a Prandti-type pitot static tube with 
a hemispherical nose behaves exactly like a Preston tube for 
boundary shear stress measurement. V.C. Patel (1965) has 
given a very accurate calibration curve for the Preston tube 
on smooth boundaries, while A.B. Hollingshead and N. 
Rajaratnam (1980) have produced calibration curves for the 
Preston tube on rough boundaries. 

In the jet tank the bed was painted, with an engraved 
grid. This grid obviously produced some uniform rougnness 
and before any shear stress measurements could be made it 
was necessary to make an estimate of the roughness. 

It has been shown (Schlichting, 1968) that the increase 
in drag due to a circular cavity is the same or less than 
that given by small protrusions. It is therefore reasonable, 
for the lack of any other information, to assume that this 
is true for grooves as well as circular depressions. 

For a uniform protuberance, a simple relationship for 


the equivalent sand roughness is; 


k. = scross-sectionelmarea A (ORL) 
Cross=SeCctiona leanea lB 


where kK,=equivalent sand roughness and areas A and B are 
defined in Figure 3.4 
For the case of a ridge of unit width and constant height, 


equation 3.1 reduces to 
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k = Aw es) 


Ww 
where Aw=breadth of ridge 
w=distance between ridges 


So, assuming that the grooves due to the grid can be 
replaced by equivalent ridges and that the grooves are 


nominally imm wide, leads to a kg value of approximately 


Ke oe Se Oe meni tn Lene natn (33) 


For a pitot static tube of external diameter of 


approximately 1.6mm, the value of 


Caley ees ibe 
kK, 0.11 a 


where d=external diameter of pitot static tube 

From the work of Hollingshead and Rajaratnam (1980), it 
can be found that the calibration curve for a d/k, value of 
15 follows Patels curve until 4p is approximately 4x10 and 
then enters a transition zone and becomes fully rough at 
Ap, approximately 1x10" 


where 4p, is defined as 


Ae = spd’ (325) 
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p=mass density, v=kinematic viscosity of the fluid 
Ap=dynamic pressure 
For the vast majority of the measurements made during 
the experimental runs, the boundary was smooth but for the 
fwouhigh velocity puns (nuns 1 land M2) earoundetnemnezZie 
the boundary was probably transitional. However because of 
the assumptions and approximations used for the analysis it 
was considered reasonable to assume that the boundary was 
smooth fror"akinmeasunemen ts: 
Using the approximation cf a smooth boundary it was 


possible to use Patel’s calibration curve for the ranges 


covered; 
Range oO eOG i ie <Fe 
——=- 1/40 0 i 
4Aov 
* * * 
ay, + 2 logy, Glob y + 4.0) (3.6) 
where x" = Log AP ue Sod hori 
nS p 10 (7) 
4ov ky 
2 2 
and Y = Lod, 9 cr tal 
0 (orc) 
Ala 
* 
Brand © | 0 <a ten 
* * ee oe 
y =a OMo2 oe) mam rs Ot emcmerra Or. 143) 7 eax - 0.0060 x 


(3.9) 


ay eas , 
| | Moan Sore < 
BOR BY 
gisiec Shs GATS at foliations Te as 


rei 7) ¢avcosd .anveras ter - 
’ ; ht gt2y (Ss ert) 107 bed : i ae 
Bew yvishraxs en) hoe rots 


Gita: aber s7nm 
Jeni «2 tua ATORMEG 


eri 7 we COLO ritoexvae ¥ ) abet 
Exons" ef? oo? Panera tes Hc ‘pase 
¢ a 


him 
cit 
= 


o 
I 


; 
: 7 7 a 
a 
¥ 
vi ao | 


; * 2 
a2) corte “yRtet) Igpeet- aos Wea 


50 


* * 
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and 4p,=Preston tube reading (the difference between 
Pitot and static pressures) 
Ty=wall shear stress 
Ah=head of water equivalent to pressure difference 
g=acceleration due to gravity 
From temperature measurements taken of the water used 
in the jet tank, it was found to be of the order of ee 
therefore the value of the Kinematic viscosity was taken as 
1.3mm°/s (Massey, 1975) 
Using the equations given above it was possible to 
construct a calibration curve, converting the readings from 
the transducer in mm of water directly into shear stress 


(N/m?}. This curve is given in Figure 3.5. 
G. Scour Hole Measurement 


A fairly elementary study of the scour hole formed by 
the jet issuing across the sandbox in the lake model was 
performed. The model was left in operation for a minimum of 
18 hours and occasionally for as long as 56 hours in order 


that the scour hole reached the asymptotic state (Rajaratnam 
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and Berry, 1977). 

Once the scour had reached the asymptotic state, the 
flume was drained and the scour hole measured along its 
central axis ,using the point gauge described in section C. 
The flume was refilled before the scour hole was filled in 
and leveled. The periodic draining of the lake model ensured 
that there was no temperature difference between the body of 
water contained in the flume and the water of the jet. This 
was also confirmed by temperature readings taken in the head 


canal and the lake model. 
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Figure 3.5 Calibration curve for bed shear stress 


INV GANAS oS OR eAPER IMENT ARSRESUETS 


APPAR Ta! j- SJE TAEXPANSION 


Introduct ion 


This section presents the results obtained from the 
photographic study made of the three dimensional turbulent 
bluff wall jet with varying finite tailwater conditions. 
Based on the present observations, a relationship for the 
change of the angle of the jet expansion is presented in a 


Suitable form. 


Review of Existing Work 


Unfortunately there appears to be no previous studies 
made on the problem of the change of expansion of a bluff 
wall jet due to variations of the tailwater depth. However 
some results have been obtained for the limiting case of an 
infinitely deep tailwater. 

Due to the problem of defining the exact position of 
the jet boundary, past investigators have only looked at the 
growth of ‘the half iwidth; ie the position wheres the loca! 
velocity is half the maximum velociy. This half width has 
been found to grow linearly with x, the distance downstream 
of the nozzle, and for wal! jets the growth is much larger 


parallel to the bed than normal to it (Rajaratnam and 
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Pani, 1970; Newman et al.,1972). Using the terminology shown 
in the definition sketch Figure 4.1 ,Newman et al. (1972) 
obtained a value for tane”’ , & being the angle of half width 
growth, of 0.278, while Pani(1972) measured a value of 0.20 
and Chandrasekhara Swamy and Bandyopadhyay (1975) arrived at 
a value of 0.166. However the experimental conditions of 
Pani were the closest to those of the present investigation 
and a value of 0.20 for tané® can be accepted. It is 
necessary to find a relationship between tan® and tane , 
where @ is the angle of growth for the nominal jet Peder’, 
Pani (1972) has shown that the velocity profiles in the 
transverse direction are similar and that the error function 
curve can be fitted to the data with reasonable accuracy, 
although the data tend to lie above the curve for larger 
Values Of zZ/by% 72 being tne; transverse dis Pancesmeacmned 
from the centre of the jet and Dey ine hakft width an sine 
transverse direction. Rajaratnam (1965), while investigating 
the modelling of the hydraulic jump by wall jets, used an 
approximate value of 2.25 for the non-dimensional length 
when the velocity became zero. Using a slightly higher 
estimate of 2.50, because the data lie just above the error 
function curve, leads to the result that the value of tané 
for an infinite tailwater is of the order of 0.50. 

The other limiting condition is that of no tailwater. 
The situation can now be described as a stream expanding 
past an abrupt expansion and as such the expansion, of the 


stream, has been found to be dependent on the Froude number 
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56 
and linear (Rouse, 1949). 
Exper iments 


It was found that the injection of dye into the jet 
produced a visual demarcation between the jet and 
surrounding fluid: From photographs taken of the ‘coloured 
jet, it was observed that the jet was basically linear in 
its expansion, and that the expansion of the nominal jet 
boundary could be measured directly and with greater 
accuracy than by velocity measurements. 

In total approximately one hundred different 
combinations of tailwater depth, initial velocity, nozzle 
aspect ratio and shape were investigated. Photographs were 
obtained and analysed for each run. In the jet tank, three 
different nozzles were used, 

1) Round nozzle, area: of 2¢54cm, aspecti#catio=i..0 

2) Rectangular nozzle, 4cm wide by 2.5cm high, aspect 
rato=1.6 

3) Small rectangular nozzle, 1.58cm wide by 0.64cm 
high, aspect ratio=2.47 

In the lake model, a study was made when the bed was 
fine sand (D,, =0.1imm) and the velocity was small enough 
that no scour was observed. A second study was made when the 
aluminum plate was placed over the sand bed. The outlet 
sizes for the two studies were 


1) Rectangular outlet, 5cm wide by 3.1icm high, aspect 
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Ratvo= 161 
2) Rectangular outlet. bcm wide by 2.00cminigm, saspect 
Gdlso=2. 50 

The initial jet velocity range covered was extensive, 
the lowest being 1.79cm/s and the largest velocity was 
265.0cm/s, although this upper velocity was for the no 
tailwater condition. The largest velocity investigated for a 
run with a tailwater was 256.1cm/s. 

The tailwater range was varied from no tailwater to a 
taiiwater to nozzle height raviove, 13.55 wheats 
approaching the infinite tailwater condition. Most of the 
tailwater to nozzle height ratios were, however, in the 
Ranger gsbeta O20: 

It should be noted that towards the downstream end of 
the plate, in the jet tank, the jet sometimes became 
unstable and oscillated. This may have been caused by the 
jet entering the deep water at the end of the tank, but while 
investigating hydraulic jumps below abrupt symmetrical 
expansions, Rajaratnam and Subramanya (1968) noticed that 
for a certain range of small tailwaters the jet tended to 
oscillate. The instability experienced in the present study 
may have similarities to that found by Rajaratnam and 
Subramanya but further study needs to be done before any 
conclusion can be reached. However due to the relatively 
short test section used in the present work any instability 


did not affect the results observed and presented here. 
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Experimental Results 


From the photographic study undertaken it was found 
Possible tTorclassify thewiet chardeleristacs imatonw ve 
separate sets. These five classifications are defined in 
Table 4.1. 

It will be noticed that Table 4.1 also contains two 
other classifications however these are the results from the 
lake model and are only separated for ease of 
identification, they are really part of Set 5. Plates 4.1 
through to Plate 4.12 show typical jet profiles from each 


Chass micats One 


Set 1 


Set 1 was the condition of no tailwater.Plate 4.1 
Gives a typical illustration of the tTlew pattern 
observed. The expansion was found to be linear and, as 
discussed previously, the rate of expansion was 


dependent on the Froude number defined as; 


Fo =Froude number 


Up=initiaisyjet velocuty 
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Table 4.1 Jet Expansion Classification 


Set Tailwater/Nozzle Comments 
height 
1 Yaa No Tailwater 
2 HA | ee 
e 0. 3kY, /Y—S 0.6 Non-monotonic 
4 ONGC Yon ee imeat pu mEaeZ, / 2 
5 OGY Ya Pinearmcdmeme=7 6/2 
6 Me /Yo~1.0 Lake model with 


metal plate 
aa Vine Voile Lake model with 
sand bed but no 


Scour 
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obtained. However from Figure 4.3 the trend was obvious. 
higher Froude numbers have a slower jet expansion, which 
also agrees with the findings of Rouse (1949) and agrees 


with general judgement. 
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Only one run was within this set and it may be an 
oddity. However from 
that that set 2 could be the 
observation appears to be a valid result and until 


proven otherwise must be treated as such. 


The results obtained from this set were very 
interesting, in that the jet did not expand in a2 
monotonic manner. In this range the tailwater was 


approximately half the nozzle depth and the top portion 


contracted again. The reason for this contraction was 
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Figure 4.3 - Variation of tanoe with Froude number for 


bluff nozzles and no tailwater 
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not obvious but may have been due to the forward 
momentum of the lower portion of the jet and possible 
side flow into the jet due to the reduced pressure in 
the high velocity region immediately before the 
contraction. Figure 4.4 shows a possible representation 
of the jet. The length of the "bulge" appeared to be 
dependent on the Froude number -,, as can be seen in 
Plates 4.3 and 4.4. In Plate 4.3 the Froude number was 
approximately one and one-half that of the jet shown in 
Plate 4.4, the position of maximum contraction does not 
appear in Plate 423 and obviously occurred tur then 
downstream. The flow pattern was obviously very complex 
and much further research would be required to 


understand the flow in this range completely. 


Set 4 


In Set 4 the study showed that the jet expanded 
linearly. However there was an abrupt initial expansion 
and therefore the line of the nominal jet boundary did 
not pass through the ends of the nozzle (ie Z#z 9/2). A 
stuay (Figure 495) Boresther variation of the mauromzata 
1/2Z, showed considerable scatter, although a trend was 
ObV1OUS@andearsStrarcqnvelime: could be, fOr d ict 


approximation, drawn through the data. 
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Set 5 


A lethepke tiwatermnto fozaietheight ratioce yi yo. 
greater or equal than one were represented by one set. 
In this set, the jet expanded linearly and the jet 
boundary passed (approximately) through the nozzle ends 
(eee Beez alee 

Photographic evidence from Newman et al.(1972) and 
Pani (1972) has shown that, for deep submergence, the 
gel haswanktexpansi one tualgonigin gS lich yenm eheomt OT 
the nozzle, this was also observed in Plate Anceen 
Unfortunately neither Pani nor Newman et al. conducted 
any analysis of their photographs nor did they provide 
any indication of tne scale for their work. From the 
present study this difference in virtual origins appears 
to be negligible and for practical purposes may be 


discounted. 
Shes der clalel i 


As noted earlier, sets 6 and 7 were covered by set 
5. Both sets were from the lake model and the tailwater 
to nozzie height ratio was unity. Set 6 gave an 
expansion angle (tané) varying from 0.14 to a maximum of 
0.235 and an average of 0.184. Set 7 gave an average 
Value fore tan vot 082.00 wi th ea variation shrom Uaioe 


0.26. The slightly larger value of tané in set / was 
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due to the additional roughness of the sand bed. 


Correlation of the Transverse Expansion 


When all the data available on the variation of tan 
with the tailwater was plotted (Figure 4.6) together in a 
dimensionless form, they were found to lie. for Y, /Y% 70.6 
,on a single curve. It was noted that the expansion angle 
(tan@) was independent of the jet velocity. because the 
results from the lake model lay well within the data 
Scatter, For relies or yy less than 0.6 the curve was 
undefined due to the non-linearity noted in the paragraph on 
set 3 and the dependency on the Froude number discussed 
previously. Figure 4.6 presents the variation of the jet 
expansion angle for different tailwater to nozzle height 
ratios and Appendix A contains the data and plots obtained 


from the photographic study. 


Summary 


From experimental observations using a dye injection 
and photographic technique, it was found that for bluff 
nozzles the jet expansion in the transverse direction was 
dependent on the taiilwater depth. It was found that the 
nozzle height could be used as a non-dimensionalising 
parameter and if used, then a single curve was obtained for 


the variation of the expansion angle. It was also found that 
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the jet could be classified in five sets, each with distinct 
Daysica | Cheracter istics: ana that tor tai lwater te mozz ic 


height ratios less than 0.6 the Froude number became another 


important parameter. 
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Introduction 


The results of velocity measurements taken along the 
centreline plane, normal to the bed and in the transverse 
plane of maximum centreline velocity are presented. From 
these observations, the velocity profiles are shown to be 
similar and relationships for the maximum velocity decay, 
growth of length scales. surface velocity and variation in 
the length of the potential core are established. When 
possible the empirical relationships are compared with 
theoretically established laws. Also the results from shear 


stress measurements are presented in suitabie forms. 


Review of Existing Work 


The earliest experimental investigations into three 
dimensional turbulent wall jets were carried out by Viets 
and Sforza (1966) and Sforza and Herbst (1970). However the 
first detailed investigations of bluff wall jets were 
performed by Newman et al. (1972) and Pani (1972), followed 
by Chandrasekhara Swamy and Bandyopadhyay (1975). The most 
complete study was performed by Pani (1972) and for this 
reason, and because his experimental arrangement was very 


similar to that of the present study, the results obtained 
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by Pani will generally be used as the asymptoptic values for 
the present relationships. A note will be made of the values 
assigned to the asymptoptic case when they occur in the 


following sections. 
Exper iments 


A total of twelve runs were made in the jet tank 
described in chapter 3. For each run, a different tailwater 
and initial jet velocity was used. The velocity measurements 
were made by using the Prandtl-type pitot static tube 
discussed in chapter 3. The shear stress measurements were 
obtained by resting the Prandtl type pitot static tube on 
the bed and using the calibration chart shown in Fig 3.5. 
The first ten runs were made with the rectangular nozzle of 
dimensions shown in Fig 4.2/‘b), while runs 11 and 12 were 
made with the 2.54cm diameter round nozzle, shown in Fig 
LPH | clit 

In runs 1 to 6 only measurements of the vertical 
centreline velocity profiles were obtained but in runs 7 to 
12, the vertical centreline velocity profile and the 
transverse velocity profile, along the plane of maximum 
centreline velocity was measured. Also for runs 7 to 12, the 
centreline shear stress and the transverse shear stress 
values were obtained. 

The tailwater was varied between 2.77cm (ie a tailwater 


LOMMOZZ) Ce aCOLM rato On mine Omdmma<imnumeot: le lcmpira 


: 
bant aeclav Giz 


@auiay af? 4 


Z ons 
ariet Jol sal me apa eotne Pry ht 
shewiisi tare th wane Agee aot .e ‘atone 
siremes.cacm vlisolev Smt feru.aew a chles Te nm 
ssi? shtesé Jovia Sout Hite? Srtt pote ic Sa! 
ataw tinenetyeejen Rows aarle yt? .£ +93 qgeris nit 
¥D ert Jefe JOTI ape teness ot! got? sat: yt ni ‘ 
> 43 of pwede, P-aro erie on a 
ie: oe “slugger efi? oi lw eben Sige efi pial ravi. 9 
ata oop (? anus of rte , fot. a of) of needa Peas 
~wotc \alxgon boue* celaneth aghd.S eet At Pa Sinan vs 
A 


isaitesw et? Ia atnesssyesen yYire 3 ofp! area My - 
oi * ora) mo) Such Derigitie eam) 2s) ew yt iooley ap ae 
ect toe sf4*oqne vy? loolav. Salliesinges isott7av oa 

mat van to ensta ant gnats .3hiiow vttoolaveay 
on: ..° of T enn Of onl’ DetueBem aaw Yilaofey 
zee iia. apd A oti, tue zeaiia serie 


78 


tailwater to nozzle depth ratio of 4.44). The initial jet 
Velocity Vapieq strom .o1.3cm/s to: 100em/s for Frunseimtem ils 
while for the runs with round nozzle, the initial jet 
velocity was of the order of 190cm/s. The significant 


details for each run are given in Table 4.2. 


Typical Vertical Velocity Profiles 


As discussed in chapter 2, it is possible to describe 
the tailwater condition as three separate types, shallow, 
moderate submergence and deep submergence. Each tailwater 
type has its own distinct characteristics. Shallow 
tailwaters have a short potential core and pass directly 
from the potential core region into the “uniform flow’ 
transition zone. Moderate tailwaters pass from the potential 
core region into a radial type decay region before entering 
the “uniform flow" transition. While deep tailwaters never 
really enter the transition into uniform flow. Figure 4.7 
(a),(b),(c) shows typical velocity profiles from each of the 
three tailwater types. In Fig 4.7 x is the distance in 
centimetres measured downstream from the nozzle, y is the 
distance measured from the bed and u is the forward 
horizontal velocity. In each profile, it can be seen that 
the maximum velocity decays the further downstream the jet 
travels and that the surface velocity approaches the maximum 


velocity for the section further downstream. Fig 4.7(a) 
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Table 4.2 Details of Experiments 
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gives a very clear illustration of the vertical velocity 
profile becoming almost uniform but still having a large 


forward velocity. 


Similarity of Vertical Velocity Profite 


It has been common among investigators of jets issuing 
into an infinite ambient fluid, to reduce the velocity 
proriles into a similar form by using the halteiength as a 
characteristic length scale. However when investigating the 
problem of a finite tailwater it was found that the 
velocity, in the vertical plane, often never reached a value 
of half the local maximum velocity. This obviously leads to 
a need to find ancther suitable scale and after a process of 
trial amd error/ it was found that the length, b, whem the 
velocity becomes 85% of the maximum velocity could be used 
Sauistactorily as the characteristic llendin. When thissscale 
was used to non-dimensionalise the vertical distance, y, and 
the maximum velocity, Cae at that section was used to 
non-dimensionalise the velocity, it was found that for each 
run the velocity profiles became similar after a downstream 
distance of approximately twelve times the nozzle height. 
Fig 4.8(a) gives a definition sketch for the scales used in 
non-dimensionalising the velocity profiles and Figure 4.9 
shows three typical reduced vertical velocity profiles. 
Furthermore when all the reduced profiles were plotted 


together on one plot, it was found that they all collapse 
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Figure 4.8. Definition sketch for scales a) vertical b) transverse 
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OnLO ai single curve. The similarity is very good Inmethe 
region close to the maximum velocity but some scatter occurs 
further away from the bed. This scatter was not unexpected 
considering the length scale used but the general trend of 
the curve, shown im Figure 45:10). 1s clear, Frome 1g74e10 (ie 
can be seen that the similarity profile lies above the curve 
of the plane turbulent wall jet on a smooth wall (the 
classical wall jet). This is due to the greater acceleration 
of the jet at distances away from the wall, discussed 
previously imvchapter 1a it Vs=interesting to mote. that for 
deeper tailwater (runs 7 and 11) the data lie closer to the 
C.W.uU. curve than the data for shallow tailwater (runs 


oodenn 0) 


Surface Velocity 


During the course of the observations, it was noticed 
that the surface of the tailwater was accelerated as the jet 
expanded into the whole depth. Plotting the ratio of the 
surface velocity to the local maximum velocity against the 
ratio of the tailwater depth to the nozzle height, shown in 
Figure?) 411, sImdtcaledmindte tor cach nunmunewe was 
considerable scatter but that generally it could be said 
thatesthe ra liouuc/ U-oadecicases Or Mliegce;n ve ratios. No 
definitive curve through the points is possible but the 
trend indicated reinforces that which intuitively should be 


expected. However, study of the velocity profiles indicated 
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Figure 4.10 - Dimensionless vertical velocity profile 


for the bluff wall jet with finite tailwater depth 
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that the surface velocity, Us, was dependent not only on the 
tailwater depth but also on the distance, x, downstream of 


the nozzle. It appeared that 


where f,=a functional relationship. 

When the above relationship was plotted, Figure 4.12 
a correlation became evident and a single curve could be 
drawn through the values obtained from all the runs. For 
values of x/Y, less than approximately 3 the curve becomes 
undefined, this is because as the tailwater decreases the 
jet expansion into the tailwater, in the vertical plane, 
becomes greater and the surface acceleration is increased 
for a region close to the nozzle. Indeed for a vee ratio 
approaching unity the expansion of the jet into the total 


tailwater depth becomes almost instantaneous. 


Velocity Decay 


The maximum velocity at any section is the relevant 
velocity scale and it is important to be able to predict the 
maximum velocity at any position. Pani (1972) has shown that 
when vA, wneresA 4s the cross-sectional area or tne moezzie, 
is used as the length scale all the maximum velocity decay 


curves for bluff nozzles collapse onto a single curve. 
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Theory also predicts that the maximum velocity is inversely 
proportional to the distance measured from the virtual 
Chiginn Ihe Virtuaiworigin 1s aftichitious ooimmteaigmiicen 
the flow velocity approaches infinity in order that the 
conservation of momentum can be maintained from a point 
source. 

In the present study when the decay of the maximum 
velocity was plotted, it was found that the length of the 
potential core varied with the tailwater depth. It was 
therefore necessary to shift the origin in the longitudinal 
direction to where the potential core vanished, ie where the 
velocity first begins to decay, Figure 4.13. When the 
velocity decay was plotted in a non-dimensional form, using 
the shifted origin, some interesting results were obtained. 
Although the opdservations exhibited some variation, it was 
found that the results for shallower tailwaters were located 
above the adjusted deeply submerged curve, thereby 
indicating that the velocity decay is less for shallow 
tailwaters. However for tailwater to nozzle depth ratios 
greater than approximately two the curve given by Pani 
(1972) is of adequately good fit. The various velocity decay 
curves are shown in Figure 4.14. 

As noted earlier, the maximum velocity was predicted by 
theory in chapter 2 to be inversely proportional to the 
longitudinal distance but that “the propor tionality ‘should 


change with the tailwater. Eq 2.38 stated 
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defining x* as the distance from the nozzle (x=0) to the 


Vir tte l-orivgin, Eq 2238 could bevrewrittenras 
u = c 


——= * 


(x-x )/ JA 


= 
= 
‘S) 


The asymptotic values given for Eq 4.3 are C=8.4 and 

x*/¥A ~2.0 (Pani,1972). When Uo/Um, VS x was plotted for each 
run (Appendix B) the curves were found to be linear. 
Correlating the calculated values for C and the distance to 
the virtual orag+n, KM wathe tmesvariattoneof taniwater in 
non-dimensional forms produced very distinct trends, Figure 
4°55 (al,ib), to which) curves could be fitted easilyeeboth 


curves differ sharply from the asymptotic value for Vays 


ratios less than about two. 
Length of potential core 


As has been stated above and in chapter 1, as the 
tailwater decreased so did the length of the potential 
core. If the length of the core is denoted by the symbol 


L, then the dimensionless length, L, may be defined as 


ie = 
a (4.4) 
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Pani (1972) found the non-dimensional length to be 
equal to 5.8. When the dimensionless lengths observed in 
the present study were compared with the variation of 
the tailwater to nozzle depth ratio (Figure 4.16),a 


single curve described the data extremely well. 


Vertical Length Scale 


The theory in chapter 2 predicted that the length 
scales in both vertical and transverse directions should 
grow linearly in their respective directions, but that the 
growth of the scales will also depend on the tailwater 
depth. In the vertical plane an inconsistency immediately 
becomes obvious because of the finite tailwater, after a 
certain distance downstream from the nozzle, according to 
theory, the vertical length scale will be greater than the 
tailwater. This obviously is an impossibility and therefore 
it should be expected that the linearity of the growth of 
the length scale will break down before this condition 
occurs. This deviation from theory occurs when the jet 
enters a stage, termed in chapter 1, as the "uniform" flow 
transition zone. 

For each run the growth of the length scale, b, was 
plotted against the longitudinal distance (in a suitable 


dimensionless form) and although the number of experimental 


i¢ 64 > v4 
alae el stat edt bed mee 


7 | 
peas 
i VW a 
M , 
¥ i; fen oprolh pig | 2 AE. oa eat r 
: i. 6 i. 7 7 
i Homie anariscerrat wanes bra tabi gat Atod ' aol 
{ io St 
tery cud enol iseaibh evi logeaass, Wier) ai, vi Ct we 
, ee eae aty: - 
ye * vives enk me treoadon. Oz) “Tt it hae bets> 22 ont aval 
{ - 4 ; 
— : + uf 18 
vad eleic He 'Té6. Ad alg fwott pe n nig ab 
"@)}gwi ' 7 avin? wit’ Ns sound awa wel 2 : cor 
. isxar ant, an7t mesntannce sing ail nil Ta 
fi: wi) ns steesy 0a! (tlw si soe hh ala ¥ sptdnay ysdo oyaert: 
B90T a7¢ brig wtf tates oT nt id 2f (avalide at Ww 
14 Alwo if. 99 xd tomaett | ont den? Bol oman ne: blue 


a a wv. 
wititags ery: igonted. hn Reed bf ts hake we ‘the 
on) a 2 see Oe i ie na For 
erat gighall evn at 


oer ee 
Lae 
pares 


oe 
a 


Shs) 


0.0 1.0 2.0 3.0 4.0 5.0 6.0 
Y./Yo 


Figure 4.16 - Growth of non-dimensional potential core 


length with increasing tailwater 


meee 
_ eH in oT 


Reade | : 


i 
; 
7 


PO et ge ie gi + fee 
{ 


7 | nity 


f. 8 6.3 Pw 38 8 
a | 


Py R102 lattneten ees 46 tiveyale 


conte 


ey, 


96 


points for each graph was somewhat limited, it was possible 
fesrepresent most plotsabyeassiragaht Wane. Pam te 72) 


gives the equation of the growth of the vertical length 


scale as 


however the definition of the length scale, b, used in the 


present study is approximately 0.588by, and therefore Eq 4.5 


becomes 


beige 5 oO S00 26) ox 


Yo Yo 


The general case of Eq 4.6, for when the tailwater becomes 


a finite depth can be represented as 


<|o>? 


where G2=gradient of the linear growth in the length scale 


b=the intersectseon, they axis, when x=0secf the 
Growtniin, theslengitheseale andi sipevlatedvioether position ot 
the wer Wwadleorig vat Tori b 
From the plots discussed above and given in Appendix B, 


it was possible to find a value of db/dx and 6/y,g for each 
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run. When these were plotted versus the tailwater to nozzle 
depth ratio, hal ars , they exhibited a very clear trend that 
was easily described by a single curve (Figure 4.17 (a) and 
(b)). It is interesting to note how the curves differ 
markedly from the asymptotic values for Vivo less than 
about two. This deviation agrees well with the observed 
increase in expansion, in the direction normal to the bed, 
of the jet discussed in chapter 1 and the shortening of the 
potential core. These observations can be explained by the 
occurrence of greater turbulent mixing near the nozzle for 
shallower tailwaters. It is also of interest that the 
virtual origin for the vertical length scale moves from 


being behind the nozzle to in front of the outlet. 


Transverse Velocity 


Typical velocity profile 


It was shown in part I that the expansion of the 
jet in the transverse plane is severely affected by the 
tailwater conditions. However there appears to be no 
reason why the velocity distribution within the jet 
should be affected and, from observations made, this was 
found to be the case. The transverse velocity profile 
was measured at a height 6 from the bed for runs 7 to 
12. 6 is the distance at which the maximum jet velocity 


occurs at the section in question, see Fig 4.8(b) for a 
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definition sketch. The velocity profiles appeared to be 
of a very similar form to those obtained from studies on 
infinite tailwaters and a typical profile obtained in 
the present observations is shown in Figure 4.18. 
Exploratory velocity readings indicated that the jet 
was, within experimental error, symmetrical and 


therefore only measurements of one half were taken. 
Dimensionless transverse velocity 


Using the half length, b , defined as in chapter 2 
and shown in the definition sketch Fig 4.8(b), as the 
length scale and using the maximum local velocity as the 
velocity scale, it was possible to reduce the transverse 
velocity profiles for each run into a single curve, 
providing the distance was approximately greater than 
eight times the nozzle height (Figure 4.19). Furthermore 
when the data points for every run were plotted 
together, it was found that they collapsed onto a single 
curve (Figure 4.20). This curve can be described quite 
well by the error function, which has been shown by 
analysis to describe the diffusion of laminar flow and 
has been found satisfactory in turbulent flow.The 
equation can be written as 
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It can be seen from Fig 4.20 that the data lie above 
the error function curve as the distance away from the 
jet centreline increases. This may be due to the 
inaccuracies caused by the velocity measuring apparatus, 
that was relatively insensitive to low velocities, as 
discussed in chapter 3. However it should be pointed out 
that this trait has also been observed by some other 
investigators (Pani,1972 ;Newman et al.,1972 ; 
Chandrasekhara Swamy and Bandyopadhyay.1$75). It may be 
that this trend is the true profile in the transverse 
direction but measurements made with a sensitive 
instrument, capable of measuring the low velocities 
occurring at the edge of the jet ,should be used to 
confirm these resuits. However the use of the error 
function curve to describe the transverse velocity 
profile is acceptable, providing it is recgonised that 
the true profile may be just above the tailend of the 
curve. It should be noted that the value of the maximum 
velocity, used for the velocity scale, is equivalent to 
Umo ana the methodsor obtaining tnesvd ile. oC seuss enor 
any given longitudinal distance, has been discussed in 


detail in section on "Velocity Decay”. 
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Pani (1972) has shown that there is considerable 
adriticultyain sobtainming a denera) equation for tthe iine 
representing the growth of the transverse length scale, 


b for ithe rease vof anmeint ini te: tanlwater epithe did 


Z y 
obtain an equation that fits the data reasonably well 


=e 2p (4.9) 


Were Zy-wid thro ie enemnozz le 


For a varying tailwater Eq 4.9 can be written in a 


more general form 


Z (4.10) 


where b, =the intersect of the z axis when x=0 

When b, versus x was plotted in a non-dimensional 
form (Appendix B}, the growth was found to be linear and 
the values for db,/dx and b,/z were calculated for 
each run. Unfortunately there were only a few values for 
both the parameters. However the change in db,/dx should 
be related to the expansion of the jet, investigated in 
Hatt wi. wUSTNOsamndcLOntOln2 .oUh. aS COlSCUSSeCHa nm iceman 
to reduce the curve given in Fig 4.5, it was found that 
the values of db,/dx obtained lay above the curve but 


within the range of values shown in Fig 4.6. When a 
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value of 2.25 was used to reduce the curve, the data lay 
about it and are shown in Figure 4.21(a}). For the 
variation of Om there was no additional information 
to enable the development of an accurate curve. 
Therefore the line shown in Figure 4.21(b) is only a 
possibility fore Viswirues jocation DUTAUne Iemore 
information becomes available it can be used as a first 


estimate. 


Centreline Bed Shear Stress 

It was shown in chapter 2 that the bed shear stress, 
Pome wdsed  rUnC GOnsOT “ne heVvno TCS numMbe haw iiomende cor 
brings added complications into any analysis and therefore 
' the results presented here can only be regarded as a 
preliminary investigation of shear stress. The measurements 
made of the shear stress are presented in Appendix B. For 
runs 7 to 10 the Reynolds number is of the same order of 
magnitude and can therefore be regarded as a constant. 
Defining the local skin friction coefficient as 


OS = tye (4.11) 


1/2uS, 


and the COetricienl Cte skin friction is 
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Some interesting observations can be made by plotting 
the coetfhicients inedassuitable sfonm.) -rom — i dumes4= 227 au. t 
can be seen that for lower tailwater runs (6,10), the 
Vena CON Sin Une OCalieCoet hic Tent Ce an Seo eas 
Significant as for the deeper runs (7,8). It is also 
interesting to note that close to the nozzle the coefficient 
is higher for the low tailwaters, indicating a greater enone 
stress and therefore greater turbulence, which agrees with 
the observation of a shorter potential core. Figure 4.22(b) 
shows a more convenient plot as it avoids the necessity of 
finding the local maximum velocity. This plot shows that the 
centreline bed shear stress decays more slowly as the 
tailwater is reduced, this again agrees with the results of 
the investigation into the velocity decay. The results 
obtained from runs 7 to 10 lie well above the curve given by 
Pani(1972) but this is due to the much larger Reynolds 
numbers used in Pani’s work. The Reynolds numbers in runs 11 
anc 12 are of Simiviar order to these of Panmihoy2)eand: it 
can be seen on Fig 4.22(b) that the coefficients of skin 
friction are in the same range as those obtained by Pani, 
however the shear stress decay is still less than the 
infinite tailwater case. 

Dimensional predictions indicate that the shear stress 


should be inversely proportional to the square of the 
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longitudinal distance. In an attempt to predict the 
variation of shear stress decay along the centreline with 
different tailwater depths, the ratio of the initial shear 
stress to the local shear stress (tie Momo ) was plotted 
against the square of the distance (x*). Although these 
relationsnips were linear, as predicted, when attempts 
(similar to those for the velocity decay} were made to 
relate the variations in slope and virtual origin to the 
change in tailwater depth no correlation was apparent. The 


results of the attempt are presented in Appendix B. 
Transverse Shear Stress 


An elementary study of the transverse shear stress was 
made but some useful results can be obtained from the 
measurements. When the shear stress profiles were plotted 
they appeared very similar to the velocity profiles 
discussed earlier. Defining the half length, be as the 
distance from the central axis to where the shear stress is 
reduced by half and using the maximum local shear stress as 
the shear stress scale, it was possible to reduce all the 
measured profiles onto one dimensionless curve. Therefore 
the shear stress profiles can be said to be similar. As 
before, with the similar transverse velocities, the error 


function defined by Eq 4.8 but now rewritten as 


Pine 
-0,.693(2Z/b 
epee 6951275 5) (4.13) 
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descibes the curve well. In this case the data lie only 
Slightly above the error function curve for distances 
removed from the central] axis. Figure 4.23(a) and (b) show a 
typical transverse shear stress profile, while Figure 4.24 


shows the collapsed similar profile. 
* 
Length scale b 


Approaching the length scale in a similar manner to 
the analysis adopted for the velocity length scale,b 
was plotted against x in a convenient non-dimensional 
form. Although the number of points for each graph was 
very limited , it could be said that the plots appeared 


linear. Writing the equation of the line as 


(4.14) 


where dx gradient 
bs intersect of line on the z axis 
It was possible to obtain estimates for the values of 
db/dx and 6*/z, for each run. Pani (1972) gave the 
approximate values for infinite tailwater as db/dx=0.09 
and 6 =0. When the data is correlated with the tailwater 


to nozzle height ratio the trends become apparent. For 
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the variation of the gradient, the correlation to the 
change in tailwater is quite good and a single curve can 
be drawn with some confidence. Unfortunately the 
Vartation Of Siiyees with the tailwater is not so well 
defined, although a trend is obvious and a curve can be 
drawn that will serve as a useful first estimate. Figure 
4.25(a) and (b) show the two relationships. As discussed 
earlier the variation of the shear stress length scale, 
tT  , Has not been successfully found. 


mo 


Summary 


From experimental observations of bluff wall jets, it 
has been found that both the vertical and transverse 
velocity profiles are similar. It was also found that the 
transverse velocity profile agreed fairly closely with the 
error function curve. The length of the potential core was 
affected by the tailwater depth and the variation was 
plotted in a convenient form. If adjustment was made for the 
potential core length, it was found that the centreline 
velocity decayed at a slower rate for shallower tailwaters. 
The maximum velocity decay was found to be inversely 
proportional to the distance measured from the virtual 
origin. The change in the velocity decay rate and the change 
in the position of the virtual origin due to the change in 
tailwater depth was presented in a convenient form. The 


length scales for both the vertical and the transverse 
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Figure 4.25 - a)Variation of transverse shear stress 
length scale gradient b)Variation of transverse 


shear stress of length scale intercept 
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velocity profiles were found to be linear and dependent on 
the tailwater. The equations of both lines were correlated 
to the tailwater depth. A preliminary investigation of the 
shear stress indicated that the centreline shear stress 
decay was reduced as the tailwater depth became smaller and 
that the transverse shear stress profile was similar. The 
shear stress length scale in the transverse direction was 
shown to be linear and that the equation of the line was 


dependent on the tailwater depth. 


belsle 69 7 w 2 


eft) to rotig 


era sissed? rae wie fed sor as 
mia gw et tron geantei raaia 9 
éou oot ioe7ttb sarevensst add at eae sees 
seu oor) sid %o nobieope ent tart brs aan t 
iq t8Iew! | st edt no 


7 ; 1 


Mee 


Grea MeliueEROSLONeSTUBLES 


Introduction 


This section presents the results of a preliminary 
study made investigating the effects on the three 
dimensional wall jet with a finite tailwater when the bed is 
allowed to scour. The shape of the scour hole is 


investigated and shown to be similar. 


Review of Existing Work 


Rajaratnam and Berry (1977) carried out studies of the 
erosion caused by circular turbulent wall jets on the beds 
of loose sand and polystyrene particles. The experiments 
were conducted in both air and water, and had an infinite 
tailwater. Rajaratnam and Diebel (1981) extended this work 
to turbulent circular wall jets with low tailwater and/or 
narrow channel width. The results obtained indicated that 
the erosion in air behaved somewhat differently than the 
erosion in water. It was also shown that the scour hole 
could be collapsed onto a single curve provided the correct 
scales were used, however the observations indicated that 
the deposition ridge needed different scales in order to 
obtain similarity. Rajaratnam and Berry (1977) defined two 
states of scour, the evolution and the asymptotic states but 


that in both states the scours were similar and that the 
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length scales were the maximum depth of scour and the 
longitudinal distance to the point of maximum scour. 
Rajaratnam and Macdougall (1981) investigated the erosion 
caused by plane turbulent water wall jets with minimum 
tailwater and were also able to show that the scour hole 
profiles was similar. In all the studies mentioned the 


particle size mean diameter was of the order of imm. 


Exper iments 


A series of nine runs were made in the lake model 
described in chapter 3. Each experiment was allowed to run 
for a considerable length of time, at least i8hrs and 
occasionally for 56hrs in order to ensure that the 
asymptotic state was reached. The tailwater depth remained 
Sonskant for-dl line runs sal ascepinivot so tcmr butane 
outlet velocity was varied from 8.88cm/s to a maximum of 
40.4cm/s. The important flow parameters are given in table 
rele 


A photographic study was made of the jet expansion and the 


centreline scour hole profile was measured. 


Results and Analysis 


Photographic observations 
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Table 4.3 Important Flow Parameters 
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From the photographs obtained it was possible to 
measure the jet expansion. Qualitively it may be said 
that when the jet enters the scour hole a great deal of 
turbulence occurs and the jet loses momentum. The jet 
then spills out over and around the region of the 
deposition of scoured material and appears, in most 
cases, to expand linearly. Plates 4.13 and 4.14 show two 
runs, one of very little scour and the other with 
considerable scour. It may be possible that the 
expansion of the jet, after leaving the scour hole, 
could be analysed as a slender wal! jet with a shallow 
tailwater but because there appears to be no information 
available on the jet expansion for slender wall jets 
with a shallow tailwater and because the observations 
presented here are limited, no attempt is made at 
predicting the expansion (refer to Appendix C for plots 


Ancuecdetue isu facia ea 


Scour hole=s tudy 


Referring to Figure 4.26 for the definition of the 
symbols used, the centreline asymptotic scour hole 
profile was measured and the results presented in Figure 
4,27 . The asymptotic or end state is virtually a steady 
state, although there are occasional bursts of 
tURDU eneecawathinh tnewscouUreiole. EhOMm 1G) 4.2 ia tb cat 
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shape. Using the maximum eroded depth, mo, as the depth 
length scale and the distance to the point of maximum 
erosion, Xo ese tnes lomavtudina! Vength scaler ar mawas 
Found that the scour hole profiles collapse onto a 
single curve, oineeee: It will be noticed that run 
Sa does iol (ies the tee Massed) cunves Very wey imee i $4 
only had a Reynolds number of 3,565, which was the 
lowest for all the runs, and it may be that the jet was 
not fully turbulent but in the transition range. As can 
De Seen Tromrht Clan oa tne eGeDOs it) Oligt Gges Delt ice ne 
scour hole did not become similar using the stated 
length scales. This was not unexpected because this fact 
had been noticed by other investigators. Due to the 
preliminary nature of the present study, no attempt was 
made to find convenient scales for the ridge. 

Following the analysis presented by Rajaratnam and 
Berry (1977) for the variation of the length scales; 


COnS 1TOCRINgG™ 4) te CoOU icaDe: Whi Glen 
ec = : [uo Vo 2 Vt??? D, Gg; “Ac; v, 8,29] (4.15) 


where Ap is the difference between the mass density of 
the bed material and that of the fluid , g is the 
acceleration due to gravity, v is the Kinematic 
viscosity, B is the width of the lake model and D is the 
mean diameter of bed material. For the present study 


D=0.1imm and the specific gravity of the sand was 2.65. 
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Using Buckinghams Pi theorem Eq 4.15 reduces to 


Ses = fy ae ; Bake aLEEA (4.16) 
Vo gDAo/p Mv D Yt Zo 
for sufficiently large Reynolds numbers , -b¥o , the 


effects of viscosity can be ignored. The tailwater depth 
is approximately equal to the outlet depth and the 

nozzle is bluff, therefore aye ayialline Yo/Y can be dropped 
from Eq.4.16. Also the width B is very large compared to 


the outlet depth and defining the densimetric Froude 
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it has been found by the other investigators that Eq 
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The longitudinal length scale, xX,,, can be predicted 
in a similar manner 
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When the two equations 4.18 and 4.19 were plotted 
it was found, Figure 4.29(a) and (b), that the points 
described reasonably well a straight line. The variation 
of can be adequately described by the line given by 
Rajaratnam and Diebel (1981) for a finite tailwater but 
unconstrained width. However the longitudinal scale has 
more scatter ane lies below the curve obtained by 
Rajaratnam and Diebel(1981), although it is possible to 


draw a straight line through the observed data. 


Summary 


It was shown that the expansion of the jet was severely 
affected by the presence of a scour hole. The scour hole 
centreline profiles were shown to be similar and that the 
length scales correlated well with the densimetric Froude 


number . 
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Figure 4.29 - Variation of the scour hole centreline 
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V. CONCLUSIONS AND RECOMMENDATIONS FOR FURTHER STUDY 


BeeGone 1usions 


Using theoretical considerations, it was concluded that 
the theory of three dimensional turbulent bluff wall jets 
with infinite tailwater depth could be applied to bluff wal! 
jets with limited tailwaters under certain circumstances. 
The theory did not apply when the tailwater was less than 
approximately one and one third times the outlet height. For 
tailwaters greater than one and one third times the outlet 
height but less than approximately four times the nozzle 
height the theory applied for a limited streamwise distance. 

From experimental studies made of the flow, it was 
found that the jet expansion was dependent on the tailwater 
depth. For very shallow or no tailwater conditions the jet 
expansion was also dependent on the Froude number but when 
the tailwater to nozzle height ratio was greater than 0.6 
the variation with the Froude number became insignificant. 
The vertical and transverse velocity profiles were similar 
in the fully developed flow and the variation of the 
respective length scales were functions of the tailwater 
depth. The surface velocity was dependent on both the 
tailwater depth and the downstream distance from the nozzle. 
It was also found that the potential core of the jet was 
shortened by reducing the tailwater, it was also observed 


that the expansion of the jet in the vertical plane was 
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increased as the tailwater depth diminished. The use of A 
to non-dimensionalise the longitudinal direction enabied the 
decay of the maximum velocity to be correlated in a suitable 
form, utilising the virtual origin. The maximum velocity 
decay was found to be reduced for lower tailwater depths. 
Shear stress measurements along the centreline indicated 
that the shear stress decayed at a slower rate as the 
tailwater was reduced. The transverse bed shear stress 
profiles were found to be similar and that the length scale 
was linear but a function of the tailwater. Elementary study 
of the bluff wall jet with erosion showed that the jet 
expansion was severely altered by the presence of a scour 
hole but that the scour holes were similar and that the 
length scales varied linearly with the densimetric Froude 


number. 


B. Recommendations for Further Study 


During the course of this study the lack of knowledge 
about finite tailwater conditions became very obvious and 
this study only answers the problems arising from a very 
limited field of exploration. To understand the flow under 
varying shallow tailwater conditions the following topics 
are recommended for future study. 

1)A detailed investigation of bluff wall jets, for all 


tai lwaters, at distances far from the nozzle. 
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2)Measurements of the transverse velocities and the 
turbulent characteristics of the flow. 

3)Detailed investigations into the cases of no tailwater and 
very low tailwater, including pressure distribution 
measurements. 

A)Investigations of the effects on the bluff wall jet due to 
bed roughness and/or a downstream channel of limited width. 
5)More experimental observations on the bed shear stress 
under varying tailwater conditions and bed roughness. 
6)Investigations into the expansion of the jet and the scour 
produced by the bluff wall jet issuing tangentially toVan 
erodable bed under varying tailwater depths and channel 
widths. 

7)A compiete detailed study of the diffusion characteristics 
of a slender wal! jet with variations of tailwater, and the 
effect of side walls on the slender jet. Especially the case 
when the tailwater is equal to the outlet depth, which can 
model the river flow into a lake. 

8)Investigations into the effect of finite tailwater 
conditions on the diffusion of the plane jet. 

9)Arising from 7, investigations of a surface jet entering a 
body of fluid with a severely limited depth. 

10) Investigations of the effect on the diffusion of the 
turbulent jet due to cross currents in the shallow 
tailwater. 

11)A11 of the above but with a jet of different temperature 


or density than that of the ambient fluid. 
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APPENDIX A - Experimental Data on the Transverse Expansion 


Figure Al 
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APPENDIX B - Velocity and Shear Stress data 
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Figure B2 Vertical length scale 
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Figure B6 Lack of correlation of shear stress decay with 
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Figure B7 Growth of length scale b* 


APPENDIX C - Jet with Scour data 
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